Ordering of the Heisenberg spin glass in four dimensions 
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Ordering of the Heisenberg spin glass in four dimensions (4D) with the nearest-neighbor Gaussian 
coupling is investigated by equilibrium Monte Carlo simulations, with particular attention to its spin 
and chiral orderings. It is found that the spin and the chirality order simultaneously with a common 
correlation-length exponent ucG = i^SG — 1-0, i.e., the absence of the spin-chirality decoupling in 
4D. Yet, the spin-glass ordered state exhibits a nontrivial phase-space structure associated with a 
continuous one-step-like replica-symmetry breaking, different in nature from that of the Ising spin 
glass or of the mean-field spin glass. Comparison is made with the ordering of the Heisenberg spin 
glass in 3D, and with that of the ID Heisenberg spin glass with a long-range power-law interaction. 
It is argued that the 4D might be close to the marginal dimension separating the spin-chirality 
decoupling/coupling regimes. 



I. INTRODUCTION 

The Heisenberg spin-glass (SG) model, or the 
Edwards- Anderson model with the isotropic Heisenberg 
exchange interaction [l[ , has been considered as the stan- 
dard model of many real SG materials In realistic 
three spatial dimensions (3D), earlier studies in the 80's 
suggested that the isotropic Heisenberg SG did not ex- 
hibit an equilibrium SG transition at any finite temper- 
ature in apparent contrast to experiments 0^ • Then, a 
proposal was made in 1992 that the model might exhibit a 
finite-temperature transition in the chiral sector, with the 
standard SG order occurring at a temperature Tsg lower 
than the chiral-glass (CG) ordering temperature Tea, 
i.e., Tsg < Tcg The occurrence of such separate 
spin and chirality transitions is now called "spin-chirality 
decoupling" Chirality is a multispin variable repre- 
senting the sense or the handedness of local noncoplanar 
spin structures induced by spin frustration. A possible 
counter-view to such a picture might be that the spin 
and the chirality order at a common finite temperature 
i3"l3|- Although there is no complete consensus, recent 
simulations point to the occurrence of the spin-chirality 
decoupling in 3D [l5l - [l9j . For example, Ref.jiot reported 
that Tsg was lower than Tcg by about 10 ~ 15%. 

To get further insight into the issue, it might be useful 
to extend the space dimension from original d = 3 to gen- 
eral d-dimensions. In d = 1, the Heisenberg SG with a 
short-range (SR) interaction exhibits only a T = tran- 
sition. In d = 2, recent calculations suggested that the 
vector SG model, either the three-component Heisenberg 
SG [2^ or the two-component XY SG [2lj, exhibited 
a T — transition but with the spin-chirality decou- 
pling, i.e., the CG correlation-length exponent vcG was 
greater than the SG correlation- length exponent vsG: 
meaning that this T = transition was characterized 
by two distinct diverging length scales, each associated 
with the chirality and with the spin. In the opposite 
limit of d 00, the model is known to reduce to the 
mean-field (MF) model, i.e., the Sherrington-Kirkpatrick 



(SK) model. The Heisenberg SK model is known to ex- 
hibit a single finite-temperature transition, with no spin- 
chirality decoupling. In high but finite d, Monte Carlo 
(MC) study by Imagawa and Kawamura suggested that 
the spin-chirality decoupling did not occur for d — 5, 
whereas the situation in d = 4 appeared somewhat more 
marginal (j^ . 

Another useful way of attacking the issue might be 
to study the one-dimensional (ID) Heisenberg SG with 
a long-range (LR) power-law interaction proportional 
to l/r"' (r the spin distance). Indeed, several studies 
both for the Ising and the Heisenberg SGs suggested 
that the ID LR SG model with a power-law expo- 
nent a might show the ordering behavior analogous to 
the d-dimensional SG model with a SR interaction [2^- 
[3H . Even a simple empirical formula relating a and d, 
d = 2/ (2(T — 1), was proposed 23], though the relation is 
only approximate. 

Recent MC calculation on the ID LR Heisenberg SG 
by Viet and Kawamura suggested that the spin-chirality 
decoupling occurred for a > cTc, but did not occur for 
a < ac, Uc being estimated numerically to be CTc — 0.8 
(29I Isoj . If one applies the approximate d — a corre- 
spondence formula quoted above [i^] , the critical dimen- 
sion below which the spin-chirality decoupling is expected 
would be dc — 3.3, suggesting that d = 4 might lie near 
the margin of, slightly on the side of the spin-chirality 
coupling regime. Of course, the above d — a correspon- 
dence formula is only approximate, and even whether 
d = 4 is greater or smaller than dc is not clear. Previous 
simulation on the 4d Heisenberg SG, which dealt with 
the linear size of L < 10, was not definitive concerning 
the occurrence of the spin-chirality decoupling in 4D 22 1 . 

Under such circumstances, the purpose of the present 
paper is first to clarify whether the spin and the chirality 
are decoupled or not in the 4d Heisenberg SG by simu- 
lating larger systems than the ones studied in Ref . [22| . 
Since d = 4 is expected to be close to the marginal dimen- 
sion concerning the spin-chirality decoupling, we wish to 
see what kind of ordering behavior is realized for the spin 
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and the chirality near the marginal dimension. 

The present paper is organized as fohows. In §11, we in- 
troduce our model and explain some of the details of our 
MC simulation. Various physical quantities are defined 
in §111. The results of our MC simulations, including the 
spin and the chiral correlation lengths, the spin and the 
chiral Binder ratios, are presented in §IV. The SG and 
CG transition temperatures are determined by carefully 
examining the size dependence of the finite-size data. In 
§V, critical properties of the spin and of the chirality are 
investigated by means of a finite-size scaling analysis. Fi- 
nally, section §VI is devoted to summary and discussion. 



II. THE MODEL AND THE METHOD 

The model we consider is the isotropic classical Heisen- 
berg model on a 4D hyper-cubic lattice, with the nearest- 
neighbor Gaussian coupling. The Hamiltonian is given 
by 



JijSi ■ Sj 

<ij> 



(1) 



where Si = {Sf,Sf, 5f ) is a three-component unit vector, 
and < ij > sum is taken over nearest-neighbor pairs on 
the lattice. The nearest-neighbor coupling Jij is assumed 
to obey the Gaussian distribution with zero mean and 
variance J^, which is taken to be unity (J = 1) in the 
following. The temperature T is measured in units of J. 

We perform equilibrium MC simulations on the model. 
The lattices are hyper-cubic lattices with N = L'^ sites 
with L = 6, 8, 10, 12, 16 and 20. We impose periodic 
boundary conditions in all four directions. Sample aver- 
age is taken over 1300, 1200, 840, 590, 430 and 256 in- 
dependent bond reahzations for L = 6, 8, 10, 12, 16 and 
20, respectively. Error bars of physical quantities are es- 
timated by sample-to-sample statistical fluctuations over 
the bond realization. 

In order to facilitate efficient thermalization, we com- 
bine the heat-bath and the over-relaxation methods with 
the temperature-exchange technique [l^l ■ For each heat- 
bath sweep we perform 11, 15, 19, 23, 31 and 55 over- 
relaxation sweeps, while the total number of temperature 
points in the temperature-exchange process are taken to 
be 35, 51, 59, 55, 55 and 60 for L = 6, 8, 10, 12, 16 
and 20, respectively. Care is taken to be sure that the 
system is fully equilibrated. Equilibration is checked by 
following the procedures of Ref. ^19]] • 



HI. PHYSICAL QUANTITIES 

In this section, we define various physical quantities 
measured in our simulations. 



For the Heisenberg spin, the local chirality at the z-th 
site and in the /z-th direction Xifi may be defined for three 
neighboring Heisenberg spins by a scalar 



Xi/j — Si+e^ • {Si X Si-e^), 



(2) 



where {fi = x, y, z, u) denotes a unit vector along the 
/i-th axis. There are in total AN local chiral variables. 

We define an "overlap" for the chirality. We prepare 
at each temperature two independent systems 1 and 2 
described by the same Hamiltonian (1) with the same 
interaction set. We simulate these two independent sys- 
tems 1 and 2 in parallel with using different spin initial 
conditions and different sequences of random numbers. 

The /c-dependent chiral overlap, q^ik), is defined as an 
overlap variable between the two replicas 1 and 2 as a 
scalar 
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i—1 fi—x.y,z,u 



(1) (2) 



(3) 



where the upper suffixes (1) and (2) denote the two repli- 
cas of the system, and fi is the position vector of the site 
i. 

The fc-dependent spin overlap, qapik), is defined by a 
tensor variable between the a and /? components of the 
Heisenberg spin, 

1 ^ 

9"Mfc) = ]vE'5ia^l?e''■'^^ («,/3 = ^,2/,^)- (4) 

i=l 

In term of the fc-dependent overlap, the CG and the 
SG order parameters are defined by the second moment 
of the overlap at a wavevector fc = 0, 



J2) „ [(|gx(0)P)] 
ICG — 



a,/3— a;,y,2 



(5) 



(6) 



where (• • • ) represents a thermal average and [•••] an 
average over the bond disorder. The CG order param- 

f2) 

eter q)jQ has been normalized here by the mean-square 
amplitude of the local chirality. 



1 ^ 

4/yE E [^^v^]' 



(7) 



i—1 fi—x,y,z,u 



which remains nonzero only when the spin has a non- 
coplanar structure locally. The local-chirality amplitude 
depends on the temperature and the lattice size only 
weakly. 

Finite-size correlation length is defined by 



1 



/ MO)')] 

2sin(fc,„/2)\/ [(g(fc,„)2)] 



(8) 
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for each case of the chirahty and the spin, ^cG and ^sG, 
where = (27r/L, 0, 0, 0) with = \km\- For the CG 
correlation length ^cGy '^6 consider two distinct defini- 
tions depending on the mutual direction between the e^- 
vector appearing in the definition of the local chirality 
(2) and the fcm-vector. When || fcm, *-e-, = we 
call the corresponding £,cG the parallel CG correlation 
length ^QQ, whereas, when ± k^, i.e., = y,z,u, we 
call the corresponding ^cG the perpendicular CG corre- 
lation length. The perpendicular CG correlation length 
^^g. is actually defined by the mean of three equivalent 
ones, each defined in the = y, z,u directions. 
The CG and the SG Binder ratios are defined by 



[(gs(0)^)] 
[(9s(0)2)]2 



(9) 



(10) 



These quantities are normalized so that, in the thermo- 
dynamic limit, they vanish in the high-temperature phase 
and gives unity in the non-degenerate ordered state. In 
the present Gaussian coupling model, the ground state 
is expected to be non-degenerate so that both gcG and 
gsG should be unity at T = 0. 



IV. MONTE CARLO RESULTS 

In this section, we present the results of our MC simu- 
lations on the 4D isotropic Heisenberg SG with the ran- 
dom Gaussian coupling. 

We show in Fig.l the temperature dependence of the 
CG and the SG correlation-length ratios, £,sg/L in (a), 
£,CgI^ in (b), and ^cq/L in (c). As can be seen from the 
figures, both the chiral ^sg/L and the spin £,cg/L curves 
cross at temperatures which are weakly L-dependent. 
Magnified views of the crossing-temperature range are 
shown in Figs.2(a)-(c) for the spin, the perpendicular chi- 
rality and the parallel chirality, respectively. 

As an other indicator of the transition, we show in 
Fig. 3 the Binder ratios for the spin (a) and for the chi- 
rality (b). The chiral Binder ratio gcG exhibits a neg- 
ative dip. The data of different L cross on the negative 
side of gcG- A magnified view of gcG in the crossing- 
temperature region is shown in Fig. 4. 

In contrast to gcG, the spin Binder ratio gsG shown 
in Fig. 3(a) exhibits no crossing in the investigated range 
of the temperature and the lattice size, monotonically 
decreasing with L. However, gsG develops a more and 
more singular shape with increasing L, a prominent peak 
appearing for larger L. 

In the L — oo limit, the Binder ratios gsG and gcG 
should satisfy here 5 — > in the high-temperature phase. 




0.3 0.35 0.4 0.45 



FIG. 1. (Color online) The temperature and size dependence 
of the spin correlation-length ratio (a), of the perpendicular 
chiral correlation-length ratio (b), and of the parallel chiral 
correlation-length ratio (c). 



and 5 = 1 at T = 0. Hence, the asymptotic form of gcG 
in the L — > cx) limit should be like the one as illustrated in 
the inset of Fig. 3(b). In fact, such a form of g is expected 
in a system with an ordered state exhibiting a continu- 
ous one-step- like replica-symmetry breaking (RSB) 331 . 
A similar form of gcG was observed in 3D [I6, 3 isj - 
Note that the one-step-like RSB discussed here is of a 
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FIG. 3. (Color online) The temperature and size dependence 
of the Binder ratio for the spin (a), and for the chirality (b). 
The inset of (b) is a behavior of gca expected in the L ^ oo 
limit. 



tion temperatures quantitatively, we plot in Fig. 5(a) the 
crossing temperatures TcrossiL) of £,sg/L versus the in- 
verse system size 1/Lav for pairs of the sizes L and sL 
with s — 2,5/3 and 5/4, where Lav = ^^L. Likewise, 

the crossing temperatures Tcross{L) of ^cg/-^^ ^cg/^ 
and gcG are plotted versus 1/Lav in Fig. 5(b). The cross- 
ing temperature Tcross{L) is expected to obey the scaling 
form, 



FIG. 2. (Color online) Magnified views of the temperature 
and size dependence of the spin correlation-length ratio (a), 
of the perpendicular chiral correlation-length ratio (b), and 
of the parallel chiral correlation- length ratio (c). 



continuous type, in contrast to the one-step RSB of a 
discontinuous type often discussed in conjunction with 
structural glasses. In the latter case, the negative dip of 
gcG should exhibit a negative divergence at the transi- 
tion, while such a negatively divergent behavior is not 
observed here. 

In order to estimate the bulk SG and CG transi- 



Tr, 



,{L; s) = Tg + CsL' 



1 

V 



(11) 



where v is the correlation-length exponent and a; is the 
leading correction-to-scaling exponent. We fit our data 
of Tcross{L; s) for the spin or for the chirality to the above 
form (11), to extract the transition temperature {Tg = 
TcG or Tsg) and the exponent 9 {9 = 9cg or 9sg) for 
the spin or the chirality. 

For the spin, we perform a joint fit of Tcross{L; s) of 
£,Sg/L for three different values of s = 2,|,|, where 
the values of Tsg and 9sg are taken common while the 
values of Cs be s-dependent. We then find an optimal fit 
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FIG. 4. (Color online) A magnified view of the chiral Binder 
ratio gcG in the crossing-temperature range. 



for TsG = 0.391(2) and 0sg 4(2) with the associated 
value, xVDOF=0.73. 

For the chirality, we have several kinds of crossing tem- 
peratures Tcross(i; s), i.e., Tcross{L\s) of ^^g/i, ical^ 
and gcG- Then, we perform a joint fit of the data 
of Tcross{J^\s) of these three kinds of Tcross{L; s), each 
with s = 2, |,|, where Tqg and 9cg are taken com- 
mon while the values of Cg be s-dependent. We then get 
TcG = 0.390(1) and Ocg = 2.4(4) with the associated 
value, xVDOF=0.51. 

One sees from these results that the spin and the chiral 
transition temperatures agree within the error bars, i,e, 
TsG = TcG within the accuracy of 1%. This observation 
strongly suggests the absence of the spin-chirality decou- 
pling in 4D, in contrast to the case of 3D where Tsg lies 
below TcG by about 10 ~ 15%. 

For the CG transition, we have another indicator, i.e., 
the negative-dip temperature Tdip{L) of the chiral Binder 
ratio gcG, which is expected to obey the scaling form, 



Tdip{L; 



T„ + cL- 



(12) 



The data of Tiiip{L) are also shown in Fig. 5(b). As can 
be seen from the figure, T^ipiL) changes its behavior with 
increasing L. It tends to decrease with L for smaller sizes, 
while it tends to increase with L for larger sizes of L > 16. 
Indeed, such a non-monotonic size-dependence of Tdip{L) 
is expected due to the following reason. For large enough 
L, the negative-dip temperature Tdip{L) should lie below 
the crossing temperature of (]gg-> ^crossiL^ • Since the 
exponents governing the asymptotic size dependence of 
Tdip{L) and Tcross{L) are 9 and \/v which satisfy the in- 
equality 9 > l/vhy definition, Tdip{L) needs to approach 
TcG from below for large enough L. Hence, a bending-up 
behavior observed in Tdip{L) for larger L is a necessary 
changeover as expected from the argument above. 
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FIG. 5. (Color online) The crossing temperatures Tcrosa 
of several quantities between the sizes L and sL are plotted 
versus the inverse system size 1/Lav with Lav = ^^L: The 
spin correlation-length ratio £,sg/L in (a); the perpendicu- 
lar chiral correlation- length ratio (,cg/L, the parallel chiral 
correlation-length ratio (J^^/L, and the chiral Binder ratio 
gcG in (b). In (b), the dip temperatures Tdip of the chiral 
Binder ratio gcG are also shown. Solid curves represent the 
fitting curves of the data on the basis of on Eq.(ll) (see the 
text for details). The spin-glass and the chiral-glass transi- 
tion temperatures are estimated to be Tsg ~ 0.391(2) and 
TcG = 0.390(1), respectively. 



Anyway, this changeover in the observed size- 
dependence of gdip{L) makes a systematic extrapolation 
of Tdip{L) difficult. Nevertheless, as can be seen from 
Fig.5(b), our data of Tdip{L) for larger L > 16 seems 
fully consistent with the Tec-value obtained above from 
the crossing temperatures. 

As mentioned above, the negative dip of gcG shown in 
Fig. 3(b) is consistent with the occurrence of a one-step- 
like RSB |ll,[ll,[ll. The corresponding spin Binder ra- 
tio gsG shown in Fig. 3(b) also develops a more and more 
singular form with a peak structure appearing for larger 
L. If one recalls the fact that gsG takes a value unity 
at T = and approaches zero above Tsg{= Tcg) in the 
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L — !> oo limit, gsc is expected to develop a negative dip 
as in the case of gcG- In the L — oo limit, this nega- 
tive dip temperature TfHp should yield TgQ. Since TgQ 
is likely to agree with Tea in 4D, a one-step-like RSB 
is expected to arise independently of the occurrence of 
the spin-chirality decoupling. In other words, in 4D, the 
Heisenberg SG is likely to exhibit a single SG transition 
without the spin-chirality decoupling. Yet, the SG (si- 
multaneously CG) ordered state is peculiar in that the 
ordered state possesses a one-RSB-like nontrivial phase- 
space structure. 



V. CRITICAL PROPERTIES 

In the previous section, we have demonstrated that, 
in 4D, the SG and the CG transitions are likely to take 
place simultaneously, i.e., Tsg = ^CG- In this section, 
we study the critical properties of the transition on the 
basis of a finite-size scaling analysis of our MC data. In 
the absence of the spin-chirality decoupling, a natural 
expectation for the critical properties is that, as usual, 
the spin is a primary order parameter of the transition. 
Then one should have vsg = vcg and risG < VCG- The 
latter corresponds to the fact that the spin is the primary 
order parameter and the chirality is the composite of the 
spin. 

We first study the critical properties of the spin by 
means of a finite-size scaling analysis of both £,sg/L and 

(2) 

QgQ. We employ the following finite-size scaling form 
with the leading correction-to-scaling term. 



L 



X{{T-TsG)L'^'''^){l + aL 



(13) 



TsG)Li/'^s«)(l + a'L-"^«), 
~ (14) 

where a and a' are numerical constants, while X and 
Y are appropriate scaling functions. The SG transition 
temperature Tsg is fixed to Tsg — 0.39 as determined in 
the previous section. 

We begin with the finite-size scaling of £,sg/T with 
i>SG and ujsG free fitting parameters. The best fit is ob- 
tained for i>sG = 1-0 and wsg^S.O. The resulting scaling 
plot is given in Fig. 6(a). Inspecting the quality of the 
plot by eyes, we put the error bars as I'sG = 1-0(1) and 
'^SG==3(1). Note that these estimates of i^sG and ljsg are 
consistent with our above estimate oidsG — '^SG + t"' — 
4(2). Next, with assuming vsG — 1 and ujsg = 3, we 

(2) 

perform the finite-size scaling analysis of q^Q to obtain 
rjsG = —0.3(1). The resulting scaling plot is shown in 
Fig.6(b). 

We also try the type of the extended finite-size scaling 
analysis proposed by Campbell et al where the scaling 



variables are chosen to take a matching between the crit- 
ical regime and the high temperature regime in order to 
get a wider scaling regime 3J]. The resulting exponent 
values turn out to be the same as those obtained above 
by the standard analysis. 
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FIG. 6. (Color online) Finite-size-scaling plots of the spin- 
glass correlation- length ratio £,sg/L (a), and of the spin-glass 
order parameter q"^Q (b), where the leading correction-to- 
scaling effect is taken into account. The spin-glass transition 
temperature is fixed to Tsg = 0.39 as determined in §1V. The 
best fit for ^sa/L is obtained with use = 1-0 and lusg = 3.0, 
while that for q^^Q is obtained with usg ~ 1-0 (fixed) and 
VSG = -0.3. 

Similar scaling analysis is also applied to the chiral 
degrees of freedom to estimate the chiral correlation- 
length exponent vcG and the chiral anomalous-dimension 
exponent rjcG- The transition temperature is fixed to 
TcG = 0.39 as determined in the previous section. The 
finite-size scaling of the chiral correlation-length ratio 
yields i^cG = 1-0(1) and ujcg = 1.7(3). We get the same 
estimates even when we use either the perpendicular or 
the transverse CG correlation-length ratio. The resulting 
scaling plot for the perpendicular one is given in Fig. 7(a). 
These estimates of vcG and ujcg are consistent with our 
above estimate of 6sg = ^SG H — — = 2.3(4). The finite- 
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size scaling of the CG order parameter q)jQ with fixing 
vcG = 1-0 and wcg = 1-7 yields r^cG = 2.4(8). The re- 
sulting scaling plot is given in Fig. 7(b). We also try the 
extended finite-size scaling analysis a la Campbell et al 
3J] . Again, as in the case of the spin, the resulting ex- 
ponent values turn out to be the same as those obtained 
above by the standard analysis. 
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FIG. 7. (Color online) Finite-size-scaling plots of the perpen- 
dicular chiral-glass correlation- length ratio (JjgI^ (8')i s-ikI of 
the chiral-glass order parameter (b), where the leading 
correction-to-scaling effect is taken into account. The chiral- 
glciss transition temperature is fixed to Tea = 0.39 as de- 
termined in §4. The best fit for i,cal^ is obtained with 



vcG = f.O, while that for q^Q is obtained with vcG 
(fixed) and rjcG = 2.4. 



f.O 



Combining the exponent estimates obtained above, we 
finally quote as our best estimates of the spin exponents, 

i^SG = 1-0 ±0.1 , r/sG = -0.3 ±0.1, (15) 

while for the chirality exponents quote 

i^CG = 1-0 ±0.1 , 77CG = 2.4±0.8. (16) 

If one applies the standard scaling or hyperscaling rela- 
tions, one can estimate other SG exponents as a ~ —2.0, 
/3sG - 0.85, 7SG - 2.3, and 6sg - 3.7, etc. 



One sees from these estimates that the correlation- 
length exponents v for the spin and for the chirality agree 
within the error bars, i.e., I'sG — ^CGi which indicates 
the existence of only one diverging length scale at the 
transition. This observation is fully consistent with the 
absence of the spin-chirality decoupling in the 4D Heisen- 
berg SG. Our data are also not incompatible with the 
relation ujsg = ^CG within the error bars. By contrast, 
the anomalous-dimension exponents satisfy the inequal- 
ity iJSG < VCG, indicating that the spin is the primary 
order parameter as usual. If one applies the scaling rela- 
tion to the CG exponents 7cg — i'^^'ilCG)i'CG, one would 
get the CG susceptibility exponent as jcG — — 0.4± 1.1. 
The estimated value of 7cg means that the CG suscep- 
tibility does not diverge, or diverges only weakly, at the 
transition. This observation is again consistent with the 
view that the primary order parameter in 4D is the spin 
and the chirality is only composite. 

The obtained CG exponents values might be compared 
with the earlier estimates by Imagawa and Kawamura on 
the same model, i.e., vsG = 1-3(2) and rjsG — —0.7(2) 
[2^ . One sees that vsg agrees with our present estimate 
within the error bars, while rjsG deviates somewhat. In 
view of the larger sizes employed in the present study as 
compared with those of ref.[23], i.e., L < 20 vs. L < 10, 
and also of larger number of independent samples, e.g., 
840 vs. 80 for L — 10, our present estimate would be 
more trustable. 



VI. SUMMARY AND DISCUSSION 

We studied equilibrium ordering properties of the 4D 
isotropic Heisenberg SG by means of an extensive MC 
simulation. By calculating various physical quantities 
including the correlation-length ratio, the Binder ratio 
and the glass order parameter up to the size as large as 
L = 20 and down to temperatures well below Tg, we 
have found that Tsg — 0.391(2) is likely to coincide with 
TcG = 0.390(1), which indicates that the spin and the 
chirality order simultaneously in the 4D Heisenberg SG, 
i.e., the absence of the spin-chirality decoupling. If Tsg 
and TcG are to differ, the distance in transition temper- 
atures should be less than 1%. We also studied the criti- 
cal properties of the transition on the basis of the finite- 
size scaling analysis. The exponents were estimated to 
be vsG = 1-0(1) and -qsc = —0.3(1) for the spin, and 
VCG = 1-0(1) and rjcG = 2.4(8) for the chirahty. Al- 
though the SG transition in 4D is usual in the sense that 
the spin is the primary order parameter, the standard 
exponent relations vsG — ^CG and rjsG < VCG being sat- 
isfied. Yet, the SG transition is somewhat unusual in the 
sense that the low-temperature SG (simultaneously CG) 
ordered state exhibits a nontrivial phase-space structure, 
i.e., a continuous one-step-like RSB. Note that the type 
of RSB is quite different from the one observed in the 
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Ising SG, or the one observed in the mean-field limit of 
both the Ising and the Heisenberg SGs. 

As mentioned in §1, possible correspondence between 
the orderings of the d-dimensional SR Heisenberg SG and 
of the ID LR Heisenberg SG with a power-law interac- 
tion has been suggested in the literature. Although this 
correspondence is by no means exact, recent numerical 
studies both on the Ising and the Heisenberg SGs sup- 
ported such correspondence. Indeed, Katzgraber et al 
proposed a formula for the d-a correspondence , a refined 
version of the one mentioned in §1 [26|, 



where r]sQ is the spin anomalous-dimension exponent of 
the d-dimensional SR system. Now, we have an estimate 
of VSG for the d = 4 Heisenberg SG as rjso — —0.3. 
Substituting this into the r.h.s. of eq.(17) and putting 
d — A, we get a — 0.79. Together with the recent 
numerical estimate of the borderline value of CTc sep- 
arating the spin-chirality coupling/decoupling regimes, 
CTc — 0.8 [IE the d-a correspondence suggests that 
the 4D lies very close to the borderline dimensionality of 
the spin-chirality coupling/ decoupling, on the coupling 
regime only slightly. Such a view on the basis of the d-a 
correspondence seems fully consistent with our present 
MC results. 

In fact, the correspondence holds also for the criti- 
cal exponents. In the d-a analogy, the exponent vsg 
of the ID LR model should be related to that of the d- 
dimensional SR model via the relation, j/5(3[lD-LR]=(i x 
:^SG [rfD-SR] [35]. Then, our 4D result suggests that the 
corresponding ID LR model should be characterized by 
the exponent vsg — 4 x 1.0 = 4. Meanwhile, ref.[3^ gave 
i^SG = 3.6(5) and i'cg = 4.0(5) for a = 0.8 so that the 
expected relation is indeed satisfied. All these results sug- 
gest that d = A probably lies fairly close to the borderline 
dimensionality of the spin-chirality decoupling/coupling, 
may even lie just at the border. 
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